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O ■ Abstract 

'w' . A pair of B°B° mesons from T(4S) decay exhibit EPR type non-local particle-antiparticle (flavor) 

correlation, ft is possible to write down Bell Inequality (in the CHSH form: S < 2) to test the non- 
(f) , locality assumption of EPR. Using semileptonic B° decays of T(4S) at Belle experiment, a clear 

violation of Bell Inequality in particle-antiparticle correlation is observed: 

S = 2.725 ± 0.167 st at ± 0.092 syst 

(N ' 
On 

1 Introduction 

m 

In Quantum Mechanics (QM), a two-particle non-separable wavefunction lead to a puzzling non-local 
correlations. This was first pointed out by by Einstein, Podolski and Rosen (EPR) 1 in 1935 to challenge 
Qh! the completeness of QM. A simplified variant was give by David Bohm using spin correlations 

£ : 1 

|^> = ^[|t>i®U>2-U)i®|T) a ] (i) 



where | 1)j (| J,)j) describes the spin state of j particle (J =1,2) with spin up (down) respectively. 
Measurement of the state on one particle, undetermined prior to the measurement, predicts with certainty 
the outcome of the same state measurement on the second particle. 

In 1964, J.S. Bell put EPR's locality principle into a testable form of an inequality (Bell Inequality) 
PJ. In 1969, Clauser, Holt, Shimony and Home 0] translate the Bell Inequality for a pair of entangled 
photons with two polarization analyzers: 

S = E{a, 0) - E(a, 0) + E(a', 0) + E{a' , 0) < 2 (2) 

where a, a' and f3, f3' are the possible angle settings on the polarizers. E(a, 0) is the correlation function 
of the polarization measurement outcome on both sides: 

E(a,0) = J d\p(X)n(a,0,X) = J d\p(\)n(a,X)n(0,\) (3) 

where fi(a, A) is the polarization measurement outcome (= ±1), determined by the local angle setting a 
and by a global hidden variable A. p(X) is a (integrable) probability distribution. Note the important 
physical assumption that the result on one side does not depend on the setting of the other side, and 
vice-versa: p(a, A) is independent of This is the locality condition in EPR. On the other hand, 
QM predicts E(a,0) = cos(a — 0). In actual experiments, due to the finite efficiency of the detectors, 
one does not have access to E{a.,0), only coincidence rates are measurable, such as i?++(a, /?) for the 
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(+1, +1) coincidence rate. Assuming that the set of detected events constitute a fair sample, one uses 
the normalized correlation function 

g) + fl__(q, (3) - (a, g) - fl_ + (q, g) 
Ma ' Pj fl ++ (a,/3)+E__(a,/3)+ J R+_(a,/3)+ J R_ + (a,^)" U 

For 5, it suffice to consider the following one parameters set (9) of settings: a = 0, a' = 29, (3 — 9 
and /3' = 30. The Bell-CHSH inequality © then becomes: 

S(9)=3E R (6)-E R (39)<2 (5) 

Fig.^, shows the QM prediction of S(8), there is a clear violation of the Bell-CHSH inequality over 
the entire range 0° < 8 < 68.5° and a maximal violation by a factor y/2 for 8 = 45°. Many experiments 
was performed, mostly in testing polarization or spin correlations 0113 El- It has never been tested with 
particle-antiparticle (flavor) correlation. 



2 Bell tests with B mesons 

The wafefunction of T(45) — > B°B° has exactly the same formalism as the pair of photons QJ: 

IV.) = ^=[\B°)i <8> \B°h - \B°)i ® |S°> 2 ] (6) 

A measurement of the flavor (B° or _B°) on one particle, undetermined prior to the measurement, will 
determine the flavor of the second particle at the same moment. Comparing to the photons, instead of 
choosing the polarizer's angle, one measure the flavor at different decay times, t a and tb, (the effect of 
mixing is the same as polarizer rotation 9 ). The main difference being that B mesons decays which 
reduces coincidence rates |Hj, e.g. 

R B °B0(t a ,t b ) = I e - 2 *'/r Be -A t /r B ^ _ e -At/r B cos(AmdA ^ . (7) 

Accordingly the correlations function is: 

E(t a , t b ) = - e -2t'/T Be -At/T B C0S (Am d At) (8) 

where Amj is the B° — B° mixing parameter, t r — t b o w t b o = 1.542ps is the mean B° decay time, 
t' = imn(t a ,tb) and At = \t a — tb\- This damping makes it impossible to violate the Bell-CHSH inequality 
(Fig. |10| . However, if one normalizes the correlation function to the undecayed pair of B° (Eq. 0J) 
then the correlation function is exactly the same as the photon case and depends only on At: 

E R (At) = cos(Am d At) (9) 

The Bell-CHSH Inequality becomes: 

S(At) = 3E R (At) - E R (3At) < 2 (10) 
Below At w I.Ttb — 2.62ps, Bell Inequality is violated, just like in photon polarization (Fig. [JJ,) 



3 Experimental Method 

This measurement is based on data sample of 78 corresponding to 80 x 10 6 T(4S*) — > BB decays 

collected at Belle detector at the KEKB asymmetric e + e~ (3.5GeV+8GeV) collider in Japan J^j- Belle 
detector consists of a silicon vertex detector (SVD), a central drift chamber (CDC), an array of aerogel 
threshold Cherenkov counter (ACC), time-of-ffight scintillator counters, a CsI(Tl) crystal electromagnetic 
calorimeter (ECL) located inside a superconducting solenoid magnet of 1.5T. An iron flux-return yoke 
outside the coil is instrumented to detect K L and to identify muons (KLM) [T3"). 

We determine the flavor of one neutral B meson by reconstructing B° — > D*~l + v and B° — > D* + l~v 
decays where denotes lepton (muon or electron). The flavor the other neutral B meson by lepton 
tagging. Because the B mesons are almost at rest and the boost of the asymmetric collider, At is 
proportional to the decay vertexes differences in z (beam) direction: At rs Az / (3^c where (3 and 7 are 
the Lorentz boost factor. 
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3.1 B° — > D* / + zv Reconstruction 

The decay chains B° -> D*-l+v,D*~ -> Z)°?t- with L>° -> JC+tt", D° -> X+tt-tt and L»° -> 
iC + 7r + 7r _ 7r (and their charge conjugate modes, same thereafter) were selected. 

All charge tracks must have at least one two-dimensional SVD hits and radial impact parameter 
\dr\ < 0.2 cm. In addition, we require that p > 0.2GeV/c for tracks in D° — > JT + 7r - 7r + channel. Charged 
kaons are identified by combining information of the TOF, ACC and dE/dx measurements in the CDC. 
7T° candidate must have mass within 0.011 GeV/c 2 of the 7r° nominal mass and momentum greater than 

0. 2 GeV/c. The photons from ir° — > 77 decays must have energy greater than 80MeV/c. 

D° candidate are selected with \MK w (Kwinr) — M D o\ < 13MeV/c 2 for K + tt~, K + ir~ir + ir~ chan- 
nels. For K+tt-t: channel, a cut of -37MeV/c 2 < (M" &(K3l ) - M D o) < 23MeV/c 2 where M D o = 
l.865GeV/c 2 was applied. Also, the 7r° with largest Dalitz weight > 10 was selected to reduce the 
combinatorial background. 

D*~ is reconstructed by combining D° and a low momentum pion. These pion track is re-fitted to 
the B° vertex to improve D* mass resolution. The cut O.UUGeV/c 2 < (M D o^ 3 -M D o) < 0.1464Ge^/c 2 
is applied. In addition, center-of-mass (CMS) momentum P£>, > 2.6GeV/c is rejected since it is beyond 
the kinematic limit for B meson decays. 

Electrons and muons are with opposite charge to the D*~ candidate. Electron identification is based 
on a combination of CDC dE/dx information, the ACC response, and energy deposit on ECL. Muons are 
identified by comparing information from the KLM to extrapolated charged particle trajectory. Their 
momentum are between 1.4 and 2.4 GeV/c in CMS. 

B° — ► D*~l + v are selected by combining D*~ candidate with muon or electron. The B° vertex 
fit must have x 2 per degree of freedom less than 50 and the fit error less than 100/zm. The angle 
in CMS between them greater than 90 degrees. In addition, they must satisfy E* B — E* u , l — \p* B \ 2 — 
\p*d*i\ 2 + 2\p b \\P£>*i\ cos (9b,d*i) — M v w where cos(#b,d*z) is the angle between p* B and p* D ,[- A cut 
cos(8b,D'i) < 1-1 was imposed. A cut on cos(6b,d»v) > 1-1 was applied where I' is the lepton with 
momentum reversed. This cut is intended for later for background subtraction |15j . 

3.2 Flavor tagging 

All other tracks not belonging the above B° — > D*~l + v selection were used for identifying the flavor of 
the second B°. A multidimensional likelihood method is used and explained in detail in ^1]. For each 
flavor decision, we assign a quantity r, which is a MC-determined flavor-tagging dilution factor from to 

1. We select events high purity (r > 0.875). Only lepton tagging sample were used to further increase 
tagging purity. The vertex fit must have \ 2 P er degree of freedom less than 50 and the fit error less than 

140^771. 

4 Background 

Three types of background have been considered: Continuum, Fake D* and wrong _D*-lepton combination. 
4.6/6~ 1 off resonance data were used for continuum background estimation. No event passed the selection. 
It can be neglected. 

Since D*~ is reconstructed from D°tt~ and D° from K + n~, K + tt~it and K + ir + ir~ir° decay chains, 
any fake D° reconstruction will contribute to the Fake D* background. We use the sideband of 

0.1560GeVy c 2 < (M D o^ g - M D o) < 0.1640GeI//c 2 

to estimate and subtract the number of events under the signal peak (Fig. 0). It is small and has small 
impact on the value of S. 

The uncorrelated D*—l background are mainly due to the reconstruction of a true D* from one B° with 
a lepton from the other B° with a smaller fraction from Charm decay leptons or fake leptons. We reverse 
the lepton momentum vector (//) in the CMS and select events where the variable cos(8b,d*i') < 1.1. 
This is to reject correlated D*l pair while selecting events with no angular correlation and is used in the 
background subtraction. 
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source 


error 


Fake D* 


0.005 


vmcor. D I 


0.030 


lepton momentum 


0.060 


mistag cut 


0.030 


Particle ID 


0.028 


Vertices Quality 


0.023 


remaining cuts 


0.042 


total 


0.092 



Table 1: Systematic Error 



After event selection and backgrounds subtraction, we are left with 3186 events. These events are sepa- 
rated into same flavor (SF) (Fig. Et-) and opposite flavor (OF) samples (Fig. [SJa). Both the correlation 
function Eu(At) and S(At) can be formed (Fig. EJ. A clear violation of Bell Inequality can be seen at 
At = (2±0.5)ps: 

5 = 2.725 ±0.167 stat , 

5 Systematic Error 

Table contains the summary of the estimated systematic error. The systematic errors for fake D* and 
uncorrelated D* I is estimated from the statistical error, with a 2a variation, of the number of events used 
in the background subtraction. The rest are estimated by varying the cut parameters by around 20% 
(larger than expected variation) from its nominal value. 

In summary, the final result for S at At = (2 ± 0.5)ps is 

5(2 ± 0.5ps) = 2.725 ± 0.167 sta t ± 0.092 syst 
which violates Bell Inequality by over 3tr. 

6 Comparison with Quantum Mechanics Predictions 

In addition to the violation of Bell Inequality, we can compare the result with the QM prediction. We 
use monte-carlo (MC) simulation to generate B°B° pair decays according to QM correlation and folded 
with detector resolutions and efficiencies. Three types of MC events were generated: B° — > D*lv (signal), 
B° — > D**lv (it has mixing, considered as signal) and B ± — > D°**lv (background). We fit these three 
MC channels to the data (Fig. HJ) which gives 3.8% of B ± background and 4.5% of B° -> D**lv. Finally, 
we plot OS, SF, En and 5 for data and for MC (Fig.[SJl. Data and MC fits remarkably well, indicating 
a good agreement between data and QM prediction. 

7 Conclusion 

Bell Inequality is a very powerful test on local hidden variable theories. It has never been tested in particle- 
antiparticle correlation using massive elementary particles. We measure the Bell-CHSH inequality using 
semileptonic B° decays. A violation greater than 3c is observed. MC simulations confirm that the result 
is consistent with QM predictions. This result will be further improved with the 140/6" 1 data sample. 
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Figure 1: Bell-CHSH Inequality S for Photon in a singlet state (a) and T(4S') -> B°B° with (b) and 
without (a) re-normalized decay probability. The shaded area violates Bell Inequality. 
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Figure 2: Fake D* background estimation: the signal and sideband region are shown. Notice the loga- 
rithmic scale. 
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Figure 3: Correlation function Er (a) and Bell Inequality S (b) after event selection and background 
subtraction. 
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Figure 4: Background level determination from MC B ± (D**) (bottom histogram), B°(D**) + B ± {D**) 
(middle) and total (top) fitted to the data (cross). 
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Figure 5: Data (crosses) and QM prediction (histogram) for Opposite Flavor (a), Same Flavor (b), Er 
(c) and S (d). 
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